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ABSTRACT 
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Context. Alfven wave dynamics in partially ionized plasmas of the solar atmosphere shows that there is indeed a cut-off wavenumber, 
i.e. the Alfven waves with wavenumbers higher than the cut-off value are evanescent. The cut-off wavenumber appears in single-fluid 
magnetohydrodynamic (MHD) approximation but it is absent in a multi-fluid approach. Up to now, an explanation for the existence 
of the cut-off wavenumber is still missing. 

Aims. The aim of this paper is to point out the reason for the appearance of a cut-off wavenumber in single-fluid MHD. 
Methods. Beginning with three-fluid equations (with electrons, protons and neutral hydrogen atoms), we performed consecutive 
approximations until we obtained the usual single-fluid description is obtained. We solved the dispersion relation of linear Alfven 
Q j waves at each step and sought the approximation responsible of the cut-off wavenumber appearance. 

I ■ Results. We have found that neglecting inertial terms significantly reduces the real part of the Alfven frequency although it never 

becomes zero. Therefore, the cut-off wavenumber does not exist at this stage. However, when the inertial terms together with the Hall 
term in the induction equation are neglected, the real part of the Alfven frequency becomes zero. 

Conclusions. The appearance of a cut-off wavenumber, when Alfven waves in partially ionized regions of the solar atmosphere are 
studied, is the result of neglecting inertial and Hall terms, therefore it has no physical origin. 
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\ 1 . Introduction the denominator of Eq. (0. Chandrasekhar d!961l l made no ex- 

rf"\ ■ ^ m , , plicit comment about this cut-off wavenumber, only pointed out 

\> . Cut-off wavenumbers m fully ionized resistive MHD have ^ when yiscQsi Qr resistivit ^ nt> Alfven waves are 

■ been reported m several classical textbooks. For instance in d d FemaQ & pi ton and Kendall & Plumpton 

P- . * e a cl ?apter °" * e ^ of the thermal ^taMity ° f a la Y er dT9§g considered the effects of finite conductivity on hydro- 
O . of fluid heated from below when a magnetic field is present, etic wayeSi and showed ^ fof k < ly we hayg timg _ 

Chandrasekhar 096B described the behavior of Alfven waves d d wayeSi whfle for k > 2y there ^ nQ waye ation 

m a viscous and resistive medium. In this case, the Alfven wave ^ ^ Furthermorei Cramer has pointed out the same ef . 

requency, u, is given y feci, showing that when the wavenumber becomes greater than 

I J i 2R„, /L, where R„, is the magnetic Reynolds number and L a ref- 

w = ± ^'V _ 4*- v_ T l) 2 ^ 1 ~ 2^ v+ ^ ' ^ erence length, the real part of the Alfven wave frequency be- 

d \ . , . . . . comes zero. Although some of the above reported textbooks 

where V A is the Alfven speed v the kinematic viscosity, r, the haye shown ^ there h a cut _ off wavenumberi none of them has 

magnetic diffusivity, and * the wavenumber. This expression investi ted the reason for its ce in resistive sing ie. flui d 

clearly points out that for a value of the wavenumber such as MHD and its physical meaning, if any. 

^ _ | ZVa On the other hand, the real part of Eq.([T} can be written as 



X 



(Or = +kr A , (3) 



V-J] 

the real part of the frequency becomes zero, and only the imag 

inary part of the frequency remains. The cut-off wavenumber with 

means that waves with a wavenumber higher than the cut-off , 

value are evanescent. Therefore, Eq. provides us with the p _ y L _ 1 , ^2^2 
cut-off wavenumber for resistive and viscous plasmas when A A V 4y2 
the single-fluid approximation is used. When viscosity is neg- 
ligible, but resistivity still exists (resistive plasma), this cut- representing a modified Alfven speed that goes to zero for the 
off wavenumber is modified by keeping only the resistivity in cut-off wavenumber, i. e. the wave ceases its propagation. 



(4) 
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Significant parts of the solar atmosphere, namely photo- 
sphere, chromosphere and prominences are only partially ion- 
ized. Collisions between different plasma species generally lead 
to the damping of waves in magnetized plasmas. Neutral atoms 
may enhance the damping of transverse waves through ion- 
neutral collision (Braginskii 1965, Khodachenko el al. 120041 
Carbonell et al. 120101 Zaqarashvili et al. 2011a,b) and may 
lead to heating of the ambient plasma (Khomenko and Collados 
120121 1. In partially ionized plasmas and in an astrophysical 
context, Kulsrud & Pearce (1969) studied the interaction be- 
tween cosmic rays and Alfven waves in the interstellar medium, 
and showed that a cut-off wavenumber appeared there as well. 
Balsara ( 1 996 1 studied MHD wave propagation in molecular 
clouds in the single-fluid approximation, reporting that Alfven 
and fast waves become completely damped when the wavenum- 
ber attains a certain value. Recently, and in the context of so- 
lar prominence oscillations, cut-off wavenumbers have been re- 
ported again. Fo rteza e t al. (120071 120081 , Barcelo et al. d20TTb 
and Soler et al. (1201 11 1 used the single-fluid approximation to 
study the damping of MHD waves produced by ion-neutral col- 
lisions in an unbounded medium with prominence physical pro- 
perties. These authors found that the cut-off wavenumber for 
Alfven waves is given by 



k=±- 



2Va cos 8 



(j] c cos 2 8 + i] sin 8) 



(5) 



with 6 the propagation angle with respect to the magnetic field, 
and T} c the Cowling diffusivity. In this case, the modified Alfven 
speed (Barcelo et al. 201 1) is given by 



T A = V AA 1- 



(jic cos 2 8 + rj sin 2 8) 2 k 2 



4V 2 cos 2 1 



(6) 



For fully ionized resistive plasmas, r\ — r\ c and, for parallel 
propagation, we recover the Kendall & Plumpton (1964} cut- 
off wavenumber. When fast and slow waves are coupled, the 
numerical value of the cut-off wavenumber for fast waves is 
similar to that of Alfven waves, but for parallel propagation 
fast waves decouple, becoming Alfven waves, and their cut-off 
wavenumber is exactly the same as for Alfven waves. Singh & 
Krishnan (2010) studied the behavior of Alfven waves in the par- 
tially ionized solar atmosphere and they also reported the cut-off 
wavenumber. 

Since time-damped transversal oscillations have been ob- 
served in the fine structure of solar filaments, Soler et al. d2009al 
I2009bb modeled this fine structure as magnetic cylinders filled 
with partially ionized plasma, which has prominence physical 
properties, and studied the time damping of transverse oscilla- 
tions produced by ion-neutral collisions and resonant absorption. 
Again, cut-off wavenumbers for fast and Alfven waves appear, 
as described by Eq.©. 

The behaviour of the cut-off wavenumber in a partially 
ionised plasma can be understood by writing its expression in 
terms of r\ and rjc- The Cowling magnetic diffusivity, t/c, is 



1,0 = h 



i i 
— h — 



a„ 



(7) 



where T en and r e , are the electron-neutral and electron-ion col- 
lisional times, respectively, g„ represents the relative density of 
neutrals, and a„ is a friction coefficient whose expression is 



2 rrii V jmij 



(8) 



with £,„ the ion-neutral collision cross-section. The Spitzer mag- 
netic diffusivity, tj, is 



An n„e 2 



1 + 1 

* P.i T ' pn 



(9) 



then, the cut-off wavenumber, k, given by Eq. (01, becomes 
2Va cos 8 



(10) 



On the other hand, as reported by Cramer (2001), the cut- 
off wavenumber, k, can be written in terms of the magnetic 
Reynolds number. In a partially ionized plasma of the solar at- 
mosphere, the Cowling magnetic diffusivity, r\c, is several or- 
ders of magnitude greater than Spitzer' s magnetic diffusivity, t], 
therefore, neglecting rj in Eq.©, we obtain 



2V A 



T]c cos 8 



(ID 



In our case, the magnetic Reynolds number, R m , could be written 
as R m = uL/t]c, where u is a characteristic velocity, which we 
could take as the Alfven speed, and L is a characteristic scale- 
length. Then, R m can be expressed as R m = VaL/t)c, and using 
Eq. (ITTb we obtain 



k = 



2R n 



Lcos 8 



(12) 



Summarizing, cut-off wavenumbers appear when the single- 
fluid approximation is used to study MHD waves in collisional 
magnetized astrophysical plasmas. In particular, this topic is rel- 
evant in connection with MHD waves in regions of the solar at- 
mosphere such as chromosphere and photosphere, or in chromo- 
spheric and coronal solar structures such as spicules or promi- 
nences. However, and up to now, an explanation for the existence 
of the cut-off wavenumber is still missing. 

Recently, Zaqarashvili et al. (1201 lal l have shown that the cut- 
off wavenumber for Alfven waves is absent in a two-fluid MHD 
approximation, where one fluid is the ion-electron gas and the 
other fluid is the gas of neutral hydrogen. Then, it is clear that 
the cut-off wavenumber appears as a result of an approximation 
and is not an intrinsic property of waves. 

Here we study the dynamics of linear Alfven waves in the 
partially ionized plasma of the solar atmosphere and seek for 
the reason for the appearance of the cut-off wavenumber. We 
start from the full multi-fluid linear equations and make consec- 
utive approximations until the usual equations for Alfven waves 
in single-fluid MHD are recovered. After each approximation, 
we derive the corresponding dispersion relations which, once 
solved, allow us to uncover the approximation responsible for 
the cut-off wavenumber appearance. 

2. Main equations 

We study partially ionized plasmas made of electrons (e), ions 
(i) and neutral (hydrogen) atoms (n). We suppose that each sort 
of species possess Maxwell velocity distribution, therefore they 
can be described as separate fluids. 

We aim to study the dynamics of Alfven waves, therefore we 
consider incompressible plasma. We also neglect the viscosity, 
the heat flux and the heat production due to collisions between 
particles. Plasma and unperturbed magnetic field are considered 
to be homogeneous. Then, linearized fluid equations for each 
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species can be split into perpendicular and parallel components 
of perturbations with regard to the unperturbed magnetic field. 
We are interested in Alfven waves, therefore only perpendicular 
components are considered. Then, the fluid and Maxwell (with- 
out displacement current) equations can be written as 



V ■ u a± = 0, 

m a n a ^^- = -e a n a \E ± + ^u aJ 



„ ^ 1 flft j. 
VxE i = — -± 
c ot 

Vx6i = —j ± , j ± 

c 



-en e (u e± — «ij_), 



(13) 
(14) 

(15) 
(16) 



where m a and n a are the mass and the number density of particles 
a, B is the unperturbed magnetic field strength, u a± and b± are 
the velocity and magnetic perturbations perpendicular to the un- 
perturbed magnetic field, E ± and j ± are the perpendicular com- 
ponents of electric field and current density, e a = +4.8 x 10~ 10 
statcoul is the charge of electrons and protons (note, that e a = 
for neutral hydrogen) and c = 2.9979 x 10 10 cm s _1 is the speed 
of light. Plasma is supposed to be quasi-neutral, which means 
n e = «;. R a± is the change of impulse of particles a due to colli- 
sions with other sort of particles, which in the case of Maxwell 
distribution in each sort of particles has the form (Braginskii 
[19651 



R a ± = - ^ a a b{u a ± - ut±), 



(17) 



where a a b — aba are coefficients of friction between particles a 
and b. 

For time scales longer than ion-electron and ion-ion collision 
times, the electron and ion gases can be considered as a single 
fluid. This significantly simplifies the equations taking into ac- 
count the smallness of electron mass with regard to the masses 
of ion and neutral atoms. Then the three-fluid description can be 
changed by a two-fluid description, where one component is the 
charged fluid (electron+protons) and the other component is the 
gas of neutral hydrogen (Zaqarashvili et al. 201 la). 

Then we may go a step further and derive the single-fluid 
MHD equations. We use the total velocity (i.e. velocity of center 
of mass) 



_ P;«ij_ + Pn«n± 

W_L — ; 

Pi + Pn 

the relative velocity 



»/J_ ~ "n_L- 



and the total density 
P - Pi + Pn, 



(18) 



(19) 



(20) 



where p a — m a n a is the density of the corresponding fluid. Then, 
Eqs. ( fOlfToT ) are rewritten as 



1 

4^ 



(V x b ± ) x B, 



dt 4-Trpgi \nen e p%£„ 



(21) 



Vxb ± - ai "+" e " w ± ,(22) 



^ = Vx(h ± xB)- ?7Vx(V x b ± ) — Vx ((V x b ± ) xB) + 

dt 4nen e 



where 



-Vxh> ± + £,Vx (w± x B) , 



: Pi/P, Zn 



n = 



Ancr Ane 2 nl 



p„lp and 



a ei + a e . 



(23) 



(24) 



The single-fluid Hall MHD equations are obtained from Eqs. 
(12 1 M23b as follows. The inertial term (the left-hand side term in 
Eq. l22l is neglected, which is a good approximation for time 
scales longer than ion-neutral collision time, but fails for the 
shorter time scales (Zaqarashvili et al. 201 la). Then w± defined 
from Eq. (l22l is substituted into Eq. ( l23l and one can obtain the 
Hall MHD equations 



P^TT = -/-(V x ft x ) x 2», 

^=Vx( Kl xB)--^ 

ot <\nen t 

+ r]cV 2 b ± , 
where 

e n B 2 



i - 



1 + 



4na;, 



(25) 

Vx((Vxi x )xB) + 
(26) 

(27) 



is the Cowling coefficient of magnetic diffusion and the second 
term in the right side of Eq. d26b is the Hall current term modified 
by electron-neutral collisions. 

The usual single-fluid MHD equations, which are widely 
used for describing of Alfven waves in partially ionized plasmas, 
are obtained from Eqs. ( 1251126b after neglecting the modified Hall 
term in Eq. 



^ = Vx(h ± x B) + j] c V 2 b ± . 
ot 



(28) 



(29) 



3. Alfven waves in partially ionized plasmas of the 
solar atmosphere 

We study the propagation of linear Alfven waves in partially 
ionized plasmas of the solar atmosphere. Here, and for the nu- 
merical calculations, we have chosen physical parameters corre- 
sponding to solar quiescent prominences (B = 10 G, n e = 10 10 
cm -3 , T = 8000 K). The unperturbed magnetic field B is di- 
rected along the z axis of cartesian frame. Next, we consider the 
Alfven wave propagation along the magnetic field, consequently, 
we perform the Fourier analysis with exp(—imt + ikz), where m 
is the wave frequency and k is the wavenumber. 

From Eqs. (12 1 M23b the following dispersion relation is ob- 
tained 

ad 2 v [1 + (1 + v)£] co - a% n [(+a + ^a>)u - 1] CO + id a 2 g n u 2 + 



+v 



±aw((- l)-a 2 <Tw 2 +£(l +£(1 +2aw) + ((a 2 - 1)£- \)co 2 ) 



P^n 



= o, 



(30) 
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where to = Tij/(kVA), t = to e jto„ a = IcVa/coi, 5 = 5 e i/to e , 
v = aj„/a e i and f = a e „/aj„. Here Va — Bj ^jAnp is the Alfven 
speed, 5 e j = a e j/{m e n e ) is the electron-ion collision frequency, 
u>i = eBlicm,) and to e = eBj(cm e ) are ion and electron giro- 
frequencies, respectively. 

The dispersion relation given by Eq. ( l30t has six different so- 
lutions. Four of the solutions are Alfven waves polarized in two 
perpendicular planes and propagating in opposite directions. The 
remaining solutions have zero real frequency at low wavenum- 
bers, so they are purely damped oscillations, but they have a low 
real frequency at high wavenumbers. Here and in the remaining 
part of the paper we consider only Alfven waves. The dispersion 
relation in Eq. ( l30t has a very involved analytical solution, there- 
fore it has been solved numerically and the solution is plotted in 
Fig. la. This plot shows that there is no cut-off wavenumber and 
it is similar to the upper plot of Fig. 2 from Zaqarashvili et al. 
(2011a). 

The next approximation is neglecting the inertial term in Eq. 
d22l . This is equivalent to the traditional single-fluid Hall MHD 
in partially ionized plasmas, and the dispersion relation is given 
by 



<5V£(i+oV-25V£ 2 [i+<r(2+<r)]w 2 - <s 4 v 2 [i+(i+v)<n 2 +^+ r A = v 



a to~+ 



+25\{l + Of„ + <5 2 v 2 (l + £'(\ - 4£) + 2£ 2 (1 + 20) 

+(5 2 v 2 (l+0 2 £ 2 +2/a(l+0^v [6v(l + (1 + v)0 + f n ] to{co 2 -\) = 

= 0. (31) 

The numerical solution of this dispersion relation is shown in 
Fig. lb. At low wavenumbers, the wave frequency is Alfvenic, 
but, when the wavenumber exceeds the value corresponding to 
the ion-neutral collision frequency {a as 0.00005 on the x axis), 
the real part of the wave frequency gradually decreases and be- 
comes very small near a ss 0.001. Since the real part of the 
frequency never becomes zero, there is no cut-off wavenumber. 
Thus, the single-fluid approach in partially ionized plasmas with 
Hall current term does not include a cut-off wavenumber. 

Eq. OTT l can be significantly simplified as £ «c 1 i.e. a en <K 
a,„ and it becomes (Pandey and Wardle 2 008t 
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to- 1 = 0, 



. Tick kVA 

V A ~ to^i 
whose analytical solution is given by 



<ikr k± WA ] _ 
y V A SitOi 



iki] c 
~Va 



kV A 

frCOi 



The real part of the frequency is 
V2 

co r = —(A - B) 
with 

A = 
and 



(32) 



(33) 



(34) 



AkWc 



\ ~ Vl + ^ + \{fcof 

y/2kv A 



4 + k 2 



B = 



Equations (l33l and d34b clearly show that the wave frequency 
always has a real part, i.e. Hall current term (the second term 
in front of to in Eq. [32l forbids the appearance of a cut-off 
wavenumber. 

Next, we neglect the modified Hall current term in Eq. (l23l . 
which gives the dispersion relation of Alfven waves commonly 
used in partially ionized plasmas: 



2 Vck 
to + i — to -1=0, 

Va 

whose analytical solution is 



-tfafc + ^4V 2 - k\ 



2V A 



(35) 



(36) 



From Eq. d36l l. and in dimensional form, the real part of the 
Alfven frequency can be written as 



m r = ±kTA 
where 



4V 2 



(37) 



(38) 



represents the modified Alfven speed for the considered case. 
From this solution, the modified Alfven speed becomes zero 
when 



k-K- 1 ^, 

Vc 



(39) 



and it determines the cut-off wavenumber for Alfven waves. 
Therefore, Alfven waves cannot propagate for wavenumbers 
higher than this cut-off value. 

The analytical solution for the real part of the frequency in 
Eq.(l36l) is plotted in Fig. lc and becomes zero near a * 0.001, 
which corresponds to the analytical cut-off wavenumber defined 
by Eq. (f39b . Furthermore, Fig. Id displays a zoom near the x 
axis of the solutions corresponding to dispersion relations given 
by Eqs. d3TT > and (l35l l. showing that only after removing the mod- 
ified Hall current term in the induction equation the real part of 
the wave frequency becomes zero. 

The condition of evanescence i. e. (k > k c ) can be written as 



k I± ^i + %„ > 2) 

to£i\tO e 8 i n j 

where 

/ 1 1 

Sin = OCin + 

\m,n,- m n n n 



(40) 



(41) 



is the mean ion-neutral collision frequency, and Eq. ( f40b can be 
approximated as 



kV A > 2<5 in f-, 

bn 



(42) 



which shows that the Alfven frequency corresponding to the cut- 
off wavenumber is lower than the mean ion-neutral collision fre- 
quency for weakly ionized plasmas (£, <k £„), and higher for 
almost fully ionized plasmas (£ » 

Hence, the cut-off wavenumber of Alfven waves in partially 
ionized plasmas appears after neglecting the inertial term (the 
left-hand side term in Eq.l22l and the modified Hall current term 
in the induction equation (Eq. 
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Fig. 1. Real part of the dimensionless wave frequency Re(a>) versus the dimensionless Alfven frequency a = kV^jiOi in partially 
ionized plasmas, where w, is the ion gyro-frequency, a) single-fluid MHD equations with inertial term (dispersion relation given by 
Eq. [30b; b) single-fluid Hall MHD equations (dispersion relation given by Eq.[3Tt: c) single-fluid MHD equations without modified 
Hall current (dispersion relation given by Eq. [33T>: d) Zoom of the solutions in panels b) (green line) and c (red line) near x axis. 



Let us estimate the Alfven frequency corresponding to the 
cut-off wavenumber for parameters typical of solar prominence 
cores (n,'=10 10 cirT 3 , n„=2xl0 10 cm -3 ). Then, the cut-off fre- 
quency exactly equals the ion-neutral collision frequency 5 !n , 
which is about 3 s . For prominence oscillations, the observed 
periods are in the range between 30 s (Balthasar et al. 1993) 
and 10-30 h (Foullon et al. 2009), therefore, since these peri- 
ods are much longer than the typical collisional timescales, the 
single-fluid approximation can be safely used when interpreting 
prominence oscillations in terms of MHD waves. 



As we stated in Sect. [TJ the above calculations are relevant 
for MHD waves in different layers and structures of the so- 
lar atmosphere. Therefore, from the expression for the modi- 
fied Alfven speed (Eq. [38] ), the numerical value for the cut-off 
wavenumber in different partially ionized layers of the solar at- 
mosphere such as the photosphere and the low and high chro- 
mosphere, can be obtained. The physical parameters for these 
regions of the solar atmosphere are taken from the FAL-1993 
model (Fontenla et al. 119931 ), and in Fig. [2] we plot the behavior 
of the modified Alfven speed, F^ versus the wavenumber, k, for 
the layers mentioned above and for quiescent prominences. This 
figure shows that in all the considered cases a cut-off wavenum- 
ber appears, whose numerical value depends, of course, on the 
different physical conditions that modify Cowling resistivity. 
Furthermore, the figure also points out that within a certain range 
of wavenumbers, the numerical value of the modified Alfven 
speed for each region or structure considered coincides (hori- 
zontal part of the plots) with the ideal Alfven speed, Va, as is 
expected from Eq. [38] Finally, we remark that the single-fluid 
approximation is valid only for waves whose wavenumbers are 
lower than the cut-off value. 




20 "-. 

q\ 1 1 ^3 lill , 

10~ 8 10~ 6 10~ 4 10~ 2 1 10 2 
k(rrf 1 ) 

Fig. 2. Modified Alfven speed, T^, versus wavenumber for the 
photosphere (dashed line), the low chromosphere (dotted line), 
the high chromosphere (dot-dash line), and quiescent promi- 
nences (solid line) 



4. Discussion and conclusion 

Consequent approximations from multi-fluid to single fluid 
MHD allow us to find the stage where the cut-off wavenum- 
ber for Alfven waves appears in partially ionized plasmas of 
the solar atmosphere. The real part of the wave frequency al- 
ways exists in full multi-fluid equations, which means that there 
is no cut-off wavenumber of Alfven waves. The real part of wave 
frequency approaches zero, but never vanishes when the inertia 
of ion-neutral relative velocity are neglected. This approxima- 
tion corresponds to the single-fluid Hall MHD and hence there 
is no cut-off wavenumber here, either. On the other hand, the 
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cut-off wavenumber appears when, after neglecting the inertial 
term, one neglects the modified Hall current term in the induc- 
tion equation. 

The appearance of Alfven waves is due to the current j x , 
which is perpendicular to b ± and k (see Eq. [161). This current 
is due to the ion inertia in an oscillating electric field E ± and it 
plays a key role in Alfven oscillations through the force j ± X B 
(Chen [1988b . In the general case, when both electrons and ions 
move freely into electric and magnetic fields, the collision of 
plasma particles and neutral atoms may lead only to damping of 
Alfven waves and waves with any frequency may propagate in 
the medium (but remember that the wave frequency should be al- 
ways lower than ion-electron collision frequency, otherwise the 
fluid approximation is not valid). When one neglects the electron 
inertia, i.e. considers the Hall MHD approximation, the elec- 
trons are frozen in the magnetic field, but ion inertia is still at 
play. In this case, the collision between ions and neutrals be- 
comes important for shorter scales and it opposes the current j ± 
at higher wavenumbers. However, oscillations again exist for all 
wavelengths due to the Hall current, but the wave frequency is 
much lower than the Alfven frequency for shorter wavelengths 
(see green lines in Fig. 1). The critical point arises when one 
neglects the Hall current in the induction equation. Then the 
collision between particles completely blocks the current j ± at 
higher wavenumbers and leads to the appearances of the cut-off 
frequency. 

Finally, there is still no cut-off wavenumber if one neglects 
the usual Hall current but retains the electron-neutral collision: 
the real part of wave frequency very closely approaches the x 
axis, but never becomes zero. This is because the electrons are 
not completely frozen into the magnetic field due to the colli- 
sion with neutrals, and therefore the ion-neutral collision cannot 
completely block the current j ± at higher wavenumbers. 

We note that a cut-off wavenumber may appear for 
Alfven waves in the stratified solar atmosphere (Murawski, 
K. & Musielak, 120101 and references therein). But this cut-off 
wavenumber has a completely different origin because it arises 
as a result of atmospheric stratification due to the solar gravity. 

In summary, the cut-off wavenumber in Alfven waves in 
single-fluid partially ionized plasmas is due to the approxi- 
mations made when proceeding from multifluid to single-fluid 
equations and is not connected to any physical process. This 
conclusion is relevant for MHD wave studies in partially ionized 
regions and structures of the solar atmosphere. 
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